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ON DOMINANCE RELATIONS AND THE STRUCTURE OF 
ANIMAL SOCIETIES: I. EFFECT OF INHERENT 
CHARACTERISTICS 


H. G. LANDAU 


COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


_ Societies are considered in which a non-transitive dominance rela- 
tion exists between every pair of members, such as the peck-right in a 
flock of hens. A one-dimensional measure of the structure of such a 
society, h, is defined, with h — 0 for equality and h = 1 for the hier- 
archy. It is assumed that each member of the society is characterized 
by an ability vector whose components depend on individual character- 
istics such as size, concentration of sex hormone, etc., but not on social 
factors such as social rank. The distribution of abilities among members 
of the society is assumed to be given by a distribution function which 
is the same for all members, and the probability that one member domi- 
nates another is given by a function of the ability vectors of the two. 

On these assumptions formulas for the expected (mean) value and 
variance of A are determined in terms of the distribution and domi- 
nance probability functions. Some special cases are calculated, especially 
that for normaity distributed abilities and dominance probability given 
by the normal probability integral. 

Several conclusions are derived. If all members are of equal ability, 
so that dominance probability is 1/2, then any sizable society is much 
more likely to be near the equality than the hierarchy; and, as the size of 
the society increases, the probability that it will be near the hierarchy 
becomes vanishingly small. If the dominance probability is a weighted 
sum of several independent components, which make up the ability vec- 
tor, then the society is less likely to be close to the hierarchy as the num- 
ber of these components increases. The hierarchy is the prevalent struc- 
ture only if unreasonably small differences in ability are decisive for 
dominance. From this it appears that the social factors, or psychological 
factors such as the previous history of dominance, which are not included 
in the present treatment, may be of great importance in explaining the 
observed prevalence of structures very close to the hierarchy in flocks of 
domestic hens. 


1. INTRODUCTION. It was observed by T. Schjelderup-Ebbe 
(1922) that between every two hens in a flock there exists a rela- 
tion known as the “peck right” which establishes the dominance of 
one hen over the other. Similar dominance relations have been ob- 
served in social groups of many other vertebrates, although not al- 
ways as clearly defined as in domestic hens and not always holding 
between every pair of members of the group. These relations have 
been intensively studied by W. C. Allee (1938, 1949) and his students 
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(Potter, 1949; Collias, 1943; Guhl and Allee, 1944) using mainly do- 
mestic hens. 

The mathematical theory of such dominance relations has been 
investigated by A. Rapoport (1949a, 1949b, 1950), and in the follow- 
ing we make use of some of the concepts introduced by him. Some 
suggestions toward a theory were also given by C. C. Lienau (1947). 

We shall be concerned with societies containing a finite number, 
n, of members with a dominance relation holding between the two 
members of every pair. A dominance relation is a binary, asym- 
metric, non-transitive relation, 7 dominates & being written j > k. 

Although the motivation and most of the applications of this 
study are from animal societies it should be remarked that there are 
other examples of dominance relations. Any tournament—chess, ten- 
nis, etc.—consisting of a single round robin with no games ending in 
draws is an example. If the tournament is a multiple round robin 
(each contestant plays more than one game with every other one), 
the member of each pair who wins the majority of the games be- 
tween them can be said to dominate the other. The baseball leagues 
are examples. Dominance relations also occur in the von Neumann- 
Morgenstern theory of games (von Neumann and Morgenstern, 
1947) where they are needed to define a solution. These writers also 
point out that the same idea, i.e., a pairwise ordering without tran- 
sitivity, occurs in the “paper-form’” in sports and races, comparisons 
of the strength of chess players in a tournament, etc. 

A complete description of the structure of a society with a domi- 
nance relation requires, of course, the statement of the n(n—1)/2 = 
(3) dominance relations between all the pairs of members of the 
society. This statement is most conveniently formulated as a matrix 
(jx) Where 

= +1 if 7 >k ) 
@=—1 ifk>j, 
Q3=0, J=1,2,--,n 


(Lienau, 1947). However, this matrix is not unique because renam- 
ing the members, that is, permuting rows and corresponding col- 
umns of the matrix, does not change the structure of the society. 
Thus the structure is given by the set of matrices which can be 
obtained from a given one by any permutation of rows and the same 
permutation of columns. 

A geometric (topological) description of the structure can also 
be given by n points with lines connecting every pair of these points 
and a direction assigned to every line. 
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Let us call the structure, when thus completely defined, the 
“dominance structure” to distinguish it from another definition of 
structure introduced below. It is easily possible to find a lower bound 


for the number of possible dominance structures. There are 2 (2) 
different matrices which can be obtained by assigning either + 1 or 
41 t0- the Cs) elements aj, which are above the principal diagonal, 
j < k; the a; below the principal diagonal are then determined. 
When the 7! permutations of n objects are applied to the rows and 
columns of any one matrix it will go into at most n!—1 of the other 


matrices. Hence there must be at least 9 (2) /n! different dominance 
structures. This number becomes very large even for moderate n. 
For 2 = 8 there are over 6,000 possible dominance structures, and 
for n = 12 there are more than 10". 


2. SCORE STRUCTURE AND HIERARCHY INDEX. Our aim is to 
characterize the structure of the society under certain assumptions 
about the properties of the members. For this purpose we do not use 
the dominance structure, but a simpler definition of structure, the 
score structure, and from it derive a one-dimensional measure of the 
structure, the hierarchy index, with which we shall be mainly con- 
cerned. 


The score structure, V, of a society is a set of nm integers 
V = (1, Vo, -*-- , Un), Where v; means that the 7th member domi- 
nates v; of the others. Any permutation of the v; does not, of course, 
alter the score structure. This definition of structure was introduced 
by Rapoport (1949a) but without the name used here. 


The score structure can be obtained from a dominance structure 
matrix by adding the +1’s in each row. A little consideration shows 
that more than one dominance structure can give the same score 
structure. The simplest example occurs for » = 5. The following 


two matrices 

Ogi F 1 ie] 
| eee a0 creel 1 i | 
eal Hoole eas | 1 and 

Sel die Seah tees ede] homed ed 

Wee ets stn O 


both have V = (3, 2, 2, 2, 1), but they cannot be equivalent under row 
and column permutation because in the first matrix 5 > 1 with 
v; == 1 and 1 dominates all the others, whereas in the second 4 > 1 
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with v, = 2, and again 1 dominates all the others. Thus V is a less 
complete definition of the structure derived from the dominance 
structure. It will not suffice for all purposes, as will be discussed in 
a later paper, but can be used for our present purpose which is to 
show the conditions under which the structure approaches the hier- 
archy. 

The hierarchy is the structure with V = (n—1, n—2,-:---, 0) 
so that the members of the society can be ordered 


1>2>3>--->n 


with each dominating all the members below it and being dominated 
py all those above. At the opposite extreme is what we call the “equal- 
ity” with 


US Vg ’ (1) 


which can occur exactly only for n odd. 

As a one-dimensional measure of the position of any society with 
respect to the extremes of equality and hierarchy we introduce the 
“hierarchy index,” 


os 12 S(%—"—=)., (2) 


1 —— a 


The hierarchy index has the range 0 to 1, the factor 12/(n? — n) 
being chosen to make h = 1 when V is the hierarchy and, of course, 
h = 0 for equality. The quantity, h, is simply a multiple of the vari- 
ance of the v; since (n — 1) /2 is the mean of the v;. 

The remainder of this paper is concerned with the investigation 
of h. We shall find that it gives a great deal of information about 
how the structure of the society depends on the properties of its mem- 
bers. Lienau (1947) thought that it would not be possible to define 
a useful one-dimensional measure of structure, but we hope to show 
that h does fulfill this function. 


3. PROPERTIES DETERMINING DOMINANCE. The actual structure 
of the society will depend on what assumptions are made about the 
properties of the members. The model used here is a generalization 
of that introduced by Rapoport (1950) and would appear to be suf- 
ficiently broad to include most of the factors making for dominance 
in animal societies, except social factors. Some excluded factors 
which limit its applicability are pointed out below. 

It is assumed that each member is characterized by an “ability 
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vector,” %; = (4j1, @jo, -+++ , Lym). These x;q measure the individual 
characteristics which make for dominance such as Size, concentyra- 
tion of male sex hormone, etc. (Collias, 1943). However, they do 
not depend on the social characteristics such as v;. 


We now introduce probability concepts. For the definitions and 
statements of the methods and theorems used we refer to H. Cramér 
(1946). It is assumed that the abilities are distributed among the 
members of the society according to the multivariate distribution 
function F(x) = F(z;),7=1,2,----,. Assuming that the dis- 
tribution is the same for all members of the society means assum- 
ing that it is homogeneous, all members being of the same breed. 
The modification needed to take into account societies of mixed 
breeds, as were used in the experiments of J. H. Potter (1949), would 
not be difficult, although the statement of the results would be more 
complicated. 


For given values of «; and x, it is assumed that the probability 
that 7 dominates k is given by a function of the two vectors p (2; , vx) 
= Pix » 1.6, 


Pr {9 > k] =p (2;5%.) =PDa, 7 4=1,2,----, 2. (3) 
Since p;, is a probability we must, of course, have 
O0<p, <1, (4) 
and since either 7 > k or k > 7, we have 
Dik + Dui = 1. (5) 


We do not attempt to discuss here the problem, or even the pos- 
sibility, of determining the functions F(x) and p;, from observations 
on societies. The latter is the identification problem which is dis- 
cussed by T. C. Koopmans and O. Reiersel (1950). It appears from 
the results below that under certain assumptions on the form of 
F(a) and pix, these functions may not be completely identifiable on 
the basis of observations on the structure of societies, that is, not 
completely determined from such observations. However, the inde- 
terminism appears to be of not a very essential nature; and more 
important, planned experiments like the paired combats carried out 
by Collias (1943) should permit direct determination of these func- 


tions. 


4. EXPECTED VALUE OF h. We denote the mean, or expected, 
value of a random variable, y, by E(y). Now 
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1a" 

DE [(v;—%)?7], (6) 
9) gaa 

where 6 = (n — 1) /2 is the mean of each v;. Since F’ (x;) is the same 

for all 7, and pj is the same function of x; and a; for all 7 and k, 

the n terms in the summation in (6) will be the same, so that 


E(h) = 


12 
ORD he j EB [ (t= %)7). (7) 


n? ay 
If @,, Wo, +++: , &n are fixed, then the probability that v, = v is 
given by 


Pr [vw=v)] => Daj, Prj, ad:010 ep ee Ds se O20 ss ’ (8) 


where 7, to 7, are any v of the n — 1 integers 2, ---- , n and jy tO In 
are the remaining n — 1 — v, and the summation is over all terms of 
this type. This is the same as the sum of the terms in the expansion 
of 


II (Daj + Djx) (9) 
j=2 


which contain v factors »,; and n — 1 — v factors p;,. The char- 
acteristic function of v,, i.e., (e+) is then 


E (e+) — JIE (pije** + Dir)» (10) 
j=2 


the characteristic function of v, — 0 is 


n-1 


— it n 
E (e#*(r»)) —=16 ( 3 ) wu (p,,;e"* ae ji) 
(11) 


and E [(v1 — %)?] is the coefficient of (it)?/2! in the expansion of 
(11) in powers of t. Now 


at at 
n 


ee (>, e? +pne 2 Jai I> (p15 + (—1)" pa) ~~] 


2*k! 


(12) 


gee a at (it)? 
=i [i+ (21; — Dj.) —— + ) +o | , 
‘ PAS I tf 
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Let 


om Dix = 9 (2; , Xx) Se US aS (13) 
Multiplying out the product in (12) we obtain 


b it * 
E (et) =1 +——_ yg, 
Riis) 
(it)? ee 
22.21 


where > indicates summation over the (n — 1) (n — 2) terms with 
j,k=2,3,----,nandj#k. Then for fixed Les 


BOGS 0) 4 fel te 9 05e) (15) 


Since the «; are actually distributed according to the distribu- 
tion function F'(x;), we have 


Pees OT RPE Sim aacr 


E (etter) 
” (16) 


= {| Hi (ne + ene *)| dF (23) dF (x2) ---- dF (@n) 


9—2 


where the integral is the Lebesque-Stieltjes integral and the integra- 
tion is over the range of the x;. Then, as in (15), 


E [(v,— 0)?] =] ie LS 059%) dF (21) dF (2-2) ---- dF (2p) 


(17) 
we (n—1) ; 3 2 
sae eee (n—2) || fon aren | ar (x,)|. 
Finally using (7) and putting 
A(a) = | g(x, uaF), (18) 
3 
E(h) -—_| IA 9) fawar (a | (19) 
n+1 


Some applications of this formula and some of the conclusions 
which can be drawn from it are discussed in Sections 6 and 7 below. 

5. VARIANCE OF h. The variance of h, o?(h), may be similarly 
obtained by use of the characteristic function. The calculation is 
more involved and the general result is rather complicated. We shall 
only outline the steps and state the final formula as follows: 
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Ys[(Se—»)"]-m 


o2(h) = E(h*) — E?(h) =| 


ie jel 
=( 12 Kes: [oy | S E[e— orm] eo 
Ni ge j,h=E 
EAT) ce . : 
(—) { n [ee — 0+ | + mm—v| 3-9] 
— F?(h). 


The quantity E[(v, — v)*] is given by the coefficient of (it)*/4! in 
the expansion of E (e*’*) in (16). This gives 


n—l1 
BC — 02 | aes [8n — 5 + 2(8n— 7) (n—2) | ae@yar(e) 
(21) 
+ (n— 2) (a) faw@arc| : 


where 
Niro =NnN(N—1)---- (n—7+1). (22) 


Similarly LE [(v1 — 0)?(v. — 6)?] is obtained from the expan- 
sion of the characteristic function of v, — 6 and v, — 0: 


ity ity 
E (e6 tse) +ita(v2-0) ) — {| i (». aii aa Oe e 2 )| 


j=2 


n 


ite ite 
| (pee? +D2e )| UF (2) OF (es) «+ dF (ttn), 


k=3 


which gives 
E[ (v1 —%)? (v2. — 6)?] 
=e { (m1)? + 2(n—2) (n? —2n —1) 
. [ A*(2) dF (a) + 2(n— 2) f f {49cm (24) 


+ 2(n—4)A(x)[A(a) + A(y)] +B(a,w)| 


| B(x, y)dF (a) dF (y) + (n= 2) co] Jawa | ‘| 
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where we have put 


B(a,y) = [ 9(@ oy, 2)aF (2). (25) 
Then from (20) we finally have 
18(n— 2) 


Be he ee 
(n+ 1) (m+ 1) ¢) 


1+ 2(n—4) | 4:(@)aF (x) 


N — 8) (2) 


2s (3n? — 15n + 20) Jaaara | ae 
2 


: | As(@aF (a) + n=) J J}soc@, nae 
(26) 
+ 2(n—4)A (x) [A (x) +AW1+ BC} 


Ble, War ear). 


This formula is too complicated to be very useful but we can 
point out two simple cases. 

If pj, = 1/2 for all 7 and k, that is, for any pair of members 
the probability of the dominance relation going in either direction 
is the same; then g;, = 0 and all the integrals in (26) vanish so that 

18 (n— 2) 
o? (h) = ————_~- 18 n’.. (27) 
(n+ 1) (n+ 1) ¢) 


In general, when the integrals do not vanish we have 
o’(h) =O(1/n), (28) 


that is, m0?(h) — constant, as n > o. 


6. APPLICATIONS. No Dominance Bias. This case, with equal 
probability of dominance in either direction for every pair, Pj. = 
4 for all j and k, occurs if the 2; are all equal, that is, when every 
member of the society has the same ability. We then have 9, = 0 
and 


3 
E(h) =——. | 29 
(h) We (29) 


10 STRUCTURE OF ANIMAL SOCIETIES 


As soon as n becomes moderately large, E(h) becomes small and 
E(h) > 0asn— o. Remembering that 0 < hf < 1, it follows, using 
the generalized Chebychev inequality (Cramér, 1946, p. 182), that 
the probability that h will differ very much from zero becomes very 
small as n increases. Among these societies equality rather than 
hierarchy would be the rule. 

This then gives for general n the behavior of societies without 
bias toward dominance, which was given in detail for n = 3, 4, 5 by 
Rapoport (1949a). We also give a table of the distribution of V 
for n = 6. This was computed by induction from Rapoport’s Table 2 


DISTRIBUTION OF V FOR n = 6, pj, = 3 


V 215. Probability of V 35h 
(Ste Secale) 720 35 
(5,4 43rd ad) 240 St 
(5, 4, 2, 2, 2, 0) 240 31 
(5,93, 3; onal, 0) 240 31 
(4, 4, 4, 2, 1, 0) 240 31 
CRY Week Ae 58) 720 oF 
(GBic 2940) 720 27 
(4 Ati ot) 80 27 
(4, 4, 3, 8, 1, 0) 720 Ck 
(5; 93 22°4> 1) 1440 23 
(4, 4, 8, 2, 2, 0) 1440 23 
(5,8 e OOF hy 1680 19 
(Aho oe lanl) 2880 19 
(23 oo OY 1680 19 
(5; Phy Pag OR Py. 2) 144 15 
(4, 745-2, 25 26k) 1680 15 
(4, 3, 3, 8, 1, 1) 1680 15 
(3) 85 8, 8,°8) 0) 144 15 
(Gn38872, 2art) 8640 11 
(433,52; 2,.2..2) 2400 7 
(35-3) 8, ub, yt) 2400 7 
ie ey 2640 3 
E(h) =8/7 


0? (h) == 12/245 
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(loc. cit.), using the fact that the probability of V = (v1, Vs .++++ Um) 
is the sum of the coefficients of terms of the form om oh ee in the 
expansion of 


LTACE ttn) 
j ,K=1 
ek 

Suppose p;, * 1/2 for some j and k, that is, there are differ- 
ences in ability making for biases toward dominance. Then it can 
be seen from equation (19) that EH(h) will be increased over its 
value in the unbiased case; a society tends to move from equality 
toward hierarchy when any bias is introduced. This was suggested 
by Rapoport, based on nm = 3. We now determine E(h) for some 
assumed forms of py, and F(x). 

Dominance Strictly Determined. At the opposite extreme from 
absence of any bias toward dominance, we might consider the case 
in which dominance is completely determined by any difference in 
ability, that is, we assume that the ability vector has a single com- 
ponent and that 

Dix —1 forz;>x,, 
(30) 
—=( -forz;<<a7,, 
and we also assume that F'(2) is continuous, so that Pr(x;=4,) = 0. 
It is clear that in this case the only possible structure for the society 
is the hierarchy, but it is interesting to see how this follows from 
the formula for E'(h). 
We have 
Gix=1 fOr v7 hes 
——] for Tage le 
then 


A(x) =f OF (2t,) ale dF (atx) 


=F (a;) —[1—F (#;)] = 2F(a;) —1, 
and 
i A*(a) dF (a) = | [4F2(2) —4F (2) +1] dF (x) = 
so from (19) 
E(h) =1; (31) 
and since h < 1, then h = 1,.0r the hierarchy, is the only possible 
structure. From (26) we can ae obtain o?(h) = 0 for this case. 


12 STRUCTURE OF ANIMAL SOCIETIES 


Linear Dominance Probability. Suppose the probability of domi- 
nance depends linearly on the ability vector. This requires that the 
range of values of the ability factors be bounded. Let 


1 (Xja) Ss (251 y Ujaz%**y Lye 
with 0 < wja < ba. Using the notation 


S(z#)=0, «<0 


Sas OSes (32) 
a Ee eo 
we put 
is Lye Wee 1 Oe : 
Dix =D WaS Ais co aee, Sli with w. 20, 5 we— 1,” (3) 
a=1 a 


so that pj, is a weighted sum of linear functions of %ja — %a With wa 
as weights. 


Then 


m 
Wa 
ITK=> — (“ja— Wee), fOrOS Via, Vee Da, 


a=1 a 


and 


bm by m Wa 
A (xj) = I rae f 2s (tia Ba) UE (aoa » +++» Sr) 
0 0 a=1 a 


Wa 
= »? kine (2ja — &a) ’ 


ad=1 a 


where &. is the mean of axa. 


| Ar(z) ar @;) = IE a jel = (2-8 | 


a=1 a 
‘ dF (xj,, wees, Xjm) 


™ /[ Wa a ™ Wa Weg oae 
=3(=.. ff tea 
b apa Dabs 
azB 


(34) 


where o,? = i i 
ere oa variance of %j2, and pag = correlation coefficient of 
“ja and “jp - 


We could have chosen the : to be . 
( ja uncorrelated, or t 
to uncorrelated variables. Agsumnlity then pag = 0 we Se ica 
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38 m Wa Ga 2 m Wa Oa 2 
E(h) =——| 1+ (n—2) 3 eo: - (35) 
70 a= Al a=1 Oa é a=1 ba 


In this case E(h) has an upper bound which is less than one. 
Since 0 < aja < ba, we have ox? < ba?/4, so that, using (33), 


if Rea Mee n+ 2 
1+ Sw? | <3 § (36) 
A a n+i1 


This result indicates that in this case also the hierarchy will not be 
the usual structure, but structures far from equality may frequently 
occur. For rectangular distributions of the xj;2, we have 


+2 
Boy <3(* Ie 
n+1 


which is not far from equality. 
If the number of factors, m, becomes large, while none of the 
weights, Wa, approaches zero, then E(h) will become small. Thus 


for equal weights 


3 
E(h) < 
wt 


MW, — Ws == = Wy, = 


Sle 


2 n+ 4m 
E(h) <— ae Oe SAS 1s Oy 1 00), 


so that if there are a large number of significant factors in the abil- 
ity vector the structure approaches the equality. 


Normal Distribution. We now suppose that the factors of the 
ability vector are each normally distributed and uncorrelated, and 
that the probability of dominance is a weighted sum of normal prob- 
ability integrals of the differences of the factors. Using the notation 


eer (ee 
G(«“z) = aie [ e 2 at 
V 2n J -~ 


for the normal probability integral, we put 


m XLja — Xa A if eK 
pa=ZanG(=—*) with wa 20, >S>wWwa=—1, (37) 
Sa 


a=1 a=1 


and 


F (a3) =1G ( ou (38) 


a=1 Ca 
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Then 


a XLja — VXxa me Lia — “La 
in=3,| 2e(**—**)_1]=-1+23m.e( =) 
Sa a=1 a 


a= 


and 


ue “Lia ae) Li LE Ura 
Ate) = f ipl —1+230.6( ——— )| mae( — 


“ja 
=—1 oS fe ; 
ce. \/ Sao eros 


using 


S b 
| G(ax + b)dG (x) =e 


, (Landau, 1950). 
4/10 


Then 


m Lia 
fawnaran= fo =a 1-43w.6 eee, 
V 8a? > ue 
m Lja m Lia 
+43 we @( +r4> WaweG | (39) 
a=1 V 8a? ae Ga" teh V/ 8a” + Caz 


2; m 
sc (eel I aa( i ). 
V Sp? + of? a=1 Ca 


Using a method similar to that for | G(az + b)dGi(s), it can be 


shown that 


ed 


= 


ie 1 1 a? 
it G? (ax) dG (x) =— + — aresin( 
oF 4 2x sl ; 


Then (39) becomes 
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| A2@ar(a;) 


=e! fa m : 1 1 ; Cae 
=1—2> wt 4> we? no ee oe 


aI a=1 JU Sas ae 2 Ca? 


+ 43>! wa Wp=— 1-23) we 


a=1 


™ 2 ye m Cae 
+> Wa | = DS 2 aresin (—“"__ ) 
a=1 


(40) 
JU a=1 Sa- 2 a? 


Pe m ; Ca m 
ee eg Arcsin |e |) using >> Ww. 1b, 


JU a=1 Sen ae 2 Ca" a=1 
and 
2(%— 2) ule Oa 
E(h) = 1 + —— 3 we aresin (—" ; (41) 
Neel fb x a=1 Sq 2 wae 


The unbiased and strictly determined cases (29) and (31) can 
be obtained as limits from (41). For the unbiased case let sa > o, 
r 


then p;, > 4, and 
° oa" 
arcsivite la), 
ee aE 2 mae 
so 


E(h) > 


am+i1 


For the strictly determined case take m = 1 and let sa > 0, then 
pix — equation (30) and 


Z ( Tae ) n 

BLOB eee a 
Scuaia Pe Cae 
soE(h) > 1. 

Here, as in the previous case, it can be seen from (41) that in- 
creasing the number of factors in the ability vector tends to reduce 
E(h). When the number of significant factors determining domi- 
nance becomes large, the structure of the society moves toward 
equality. 

7. SIGNIFICANCE FOR FLOCKS OF HENS. The observations and 
experiments on flocks of hens by Schjelderup-Ebbe and by Allee and 
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coworkers all show that stable flocks of hens almost always have a 
structure that departs very little from the hierarchy. A flock of ten 
to twenty hens will normally have not more than two or three cycles, 
ie, 7 >k,k>landl >j. This means that h is normally very close 
to unity in sucha flock. A typical example quoted by Schjelderup-Ebbe 
is the following score structure for a flock of ten hens: V = (8, 8, 8, 
6, 5, 4, 3, 2,1, 0). This gives h = .975. 

We wish to show that if the probability of dominance depended 
only on inherent individual characteristics as assumed in the pres- 
ent treatment (and not on social factors), then the occurrence of 
societies with h near one would be unusual rather than the rule. 

We use the results for the normal distribution (Sec. 6) as be- 
ing probably closest to reality. As noted in Section 6, if the ability 
vector contains several uncorrelated factors this will tend to reduce 
the expected value of h, so that we can consider the case of a single 
factor. Then 


E(h) = 


+ Rey 
Tonia Weer Pro. 


where 


; V2 c- 
FIGURE Ly Normal distribution. Dependence of probability of dominance 
Pj,» On difference in ability. Values on curves give h,. 
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6 ; o 
ha =— arcsin ————. (42) 
7 Papa oe 


For even moderate n, E'(h) is close to h,, its asymptotic value, 
and in any case 


3 
E(h) —h,=—— (1—E(h)), 
n—2 


so we consider how h, depends on s/c. 
In Figure 1 we have drawn the probability of dominance, 


4 me 
Pa=G(— ), 


§ 


as a function of the difference in ability for various values of s/c 
with the corresponding values of h, noted on each curve. The ab- 
scissa is 

Lj — Ue 

V2e 

since the variance of x; — x, is 2 o?. It can be seen that for h, to be 
very close to one, very small differences in ability would have to be 
quite decisive as to dominance. 

This very close dependence of dominance on ability hardly seems 
reasonable, but there is also some experimental evidence on this 
point. Collias (1943) staged 200 combats between hens in which he 
measured degree of moult, comb size, weight and rank in own flock, 
and obtained for the correlation with success the values: .580, .593, 
A474, and .262. The correlation measure used was Pearson’s coefficient 
of biserial correlation, 7 , (Pearson, 1909) except for degree of moult 
which was not measured on a continuous scale. In the present case 


“1 Wy — Lo. 
Pe ee (43) 
2 « 


21-02 


where x; — #2 is the mean of x, — 2 when % > %, and o,_,, is the 
standard deviation of x, — #2. It is not difficult to determine the 
value of r in terms of s/o for our assumed normal distribution. The 


result is 


g2 yA 
ra(1t+s—) (44) 
oO 


and this gives 
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6 C 
h,=— aresin—. (45) 
cr Zz 


Hence the largest value of r obtained by Collias, 593, gives only 
h, = .34. Figure 2 shows the dependence of h, on 7. It is apparent 
that h, close to one requires r to be unreasonably large. 


6 ewe ee 
Fp aresin = 


ha= 


0 2 4 6 23 1.0 


FIGURE 2. Dependence of h, on biserial correlation coefticent, 7. 


We must conclude that factors omitted from the present treat- 
ment must be included to account for the observed frequency of 
structures near the hierarchy in flocks of hens. The most obvious 
omissions are the social factors, such as social lag, or the effect of 
existing differences in social rank on the probability of continued 
dominance. All observers agree that such social factors are of great 
importance. An attempt will be made to treat them mathematically 
in a later communication. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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A NOTE ON THE McCULLOCH-PITTS NEURAL NET FOR 
HEAT-COLD DISCRIMINATION 


HERMAN T. EPSTEIN AND ANATOL RAPOPORT 
THE UNIVERSITY OF PITTSBURGH} AND COMMITTEE ON MATHEMATICAL BIOLoGy, 
THE UNIVERSITY OF CHICAGO 


An error is pointed out in the neural net model which had been 
previously proposed to explain the heat-cold responses that depended on 
the duration of the stimulus. A corrected net for the effect is given. 


In their paper on the Boolian algebra of neural nets, W. S. Mc- 
Culloch and W. Pitts (1943) gave as an illustration a neural net 
which was to satisfy the following physiological situation: 


a) When a cold object momentarily touches the skin, a sensa- 
tion of heat develops; 

b) When the cold object is held on the skin for a longer time, 
the sensation is one of cold, with no accompanying sensation of heat. 


The authors derived the circuit of Figure 1 to represent the 
above situation. 


The simultaneous firing of two end-feet is required to fire the 
neuron on which they synapse. That is, for example, the neuron N, 
will fire provided N, fires one time interval previously and simultane- 
ously No fires one time interval previously. In the notation of N. 
Rashevsky (1948) this is expressed as 


Na(t) = Na(t—-1) -N2(t—1). (1) 


2 ee 
eae ae 


FIGURE 1 


+Contribution no. 8 of the Department of Biophysics of the University of 
Pittsburgh. ' 
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The looping connection of Nz on N; is one of complete inhibition. 
If N.2 fires at any time, then N, will not fire at one time interval 
later, no matter what its other stimulations may be. The circuit of 
Figure 1 was derived to fit the equations, 


N,(t) = No (t—1) - N2(t—2) (2) 
N;(t) = N,(t—1) VIN2 (t—3) - N2(t—2)], (3) 


where N, is the receptor for heat, N2 the receptor for cold, Ns the 
neuron which leads to centers giving a sensation of heat, and Nz 
correspondingly leads to a sensation of cold. 

The purpose of the present note is to point out that the circuit 
does not satisfy the condition that there be no sensation of heat if 
the cold object is held on the skin for more than one time interval. 
This is so, because when N, stops firing, say at (t + k), N, still fires 
at (t + k), and, since N, is now not inhibited by N., N, fires at 
(¢+k+1) and N, (the heat neuron) fires at (¢ + k + 2). 


FIGURE 2 


The net in Figure 2 is proposed as realizing conditions a) and 
b) above where b) is strengthened to 

b') When the cold object is held on the skin for a longer time, 
the sensation is one of cold, with no sensation of heat either preced- 
ing or following. 


The inhibitory connections from N; and N ; prevent the firing of 


N; by Nz at all times, except when Nz fires exactly once, as can be 
readily verified. 
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A NOTE ON THE EQUATIONS OF CONDITIONED REFLEX 


JOHN Z. HEARON 
DEPARTMENT OF PHYSIOLOGY 
THE UNIVERSITY OF CHICAGO 


An alternative method is suggested for integrating a certain dif- 
ferential equation associated with a conditioning process, where the 
stimulus is presented in the torm of a “square wave,” i.e., is of constant 
intensity during an interval of time followed by no stimulus during the 
next interval, etc. 


A solution is also given where the stimulus is a rectified sine wave. 


In discussing a simple, illustrative case of conditioning based 
on the two-factor theory N. Rashevsky considers that the conditioned 
and unconditioned stimuli are presented together at regular inter- 
vals. The excitatory factor, e, then obeys the equation 


de (t) 
— = AE (t) —ae(#) (1) 
dt 


wherein the total intensity of excitation, Z(t) has, beginning at t=—0, 
the value E for an interval ¢t,, is then zero for an interval ¢, , then 
has the value # for an interval f,, etc. 

The solution of equation (1) under the above conditions on E(t) 
was obtained by Rashevsky essentially as follows: The integral of 
(1) with E(t) = E describes e(t) during a ¢,-interval. During a 
t,-interval e(t) is described by the integral of the equation 


de (t) 
dt 


which gives the decay of « when E(t) = 0. For any period the 
initial conditions are given by the terminal conditions of the pre- 
ceding period. By piece-wise integration in this manner the condi- 
tions for an induction are set up: It is shown that if the formula 
obtained for e at the end of the second ¢,-period is true when gen- 
eralized for the nth t,-period it is then true for the (n + 1)st t,-pe- 
riod. The general formula for « at the end of the xth ¢,-period is 
thereby established. 


—=— ae(t) (2) 
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The biophysical limitations of the simple theory of condition- 
ing which gave rise to the above problem as well as the modifica- 
tions and extensions which make it more realistic have been dis- 
cussed by N. Rashevsky (1948, p. 436-42). 


It is the purpose of this note to discuss the purely mathematical 
point that the straightforward integration of (1), subject to the 
above described conditions on E(t), leads at once to the general re- 
sult given by Rashevsky. The increase in e when E(t) = E and the 
decay of « when E(t) = 0 are automatically accounted for in the 
method here described. 


The solution of (1), when «(0) = 0, is 
t t 
6) Se { eur B(A)d— A { oun =e (3) 


For the purpose of illustration consider E(t) to be the half-wave 
rectification of the square wave of amplitude E and period 2c. Of 
the two quadratures indicated in (3) the second is simpler to per- 
form; the first more clearly exhibits the features of the general meth- 


od. The first quadrature in (3) gives for « at the end of nm half pe- 
riods: 


é (2k+1)e 
e(me) = AES i e-aine-dy (4) 
k=0 2ke 


where €, the upper index of the summation is 


§ = (n—1) /2 when nis odd, 
§&= (n—2) /2 when nis even. 


When the integrations indicated in (4) are carried out the result is 


AE & 
e(nc) =—— (1—e*) § ear rake, (5) 

a k=0 
When & takes on the values 0, 1, 2 , ‘s+: (n—1)/2 the quantity 
(n—1—2k) takes on the values (n—1), (n—3), (n—5), ----, 0. When 
k takes on the values 0,1,2,----, (n—2) /2 the quantity (n—1—2k) 


takes on the values (n—1), (n—3), (n—5), «1. A sj 
Srila aMetia | ), +++», 1. A simple change 


AE (n-1)/2 
e(nc) = a (1 — e**) > pees. “arodd (6) 
1=0 


and 
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A EH (n-2) /2 
e(nc) =—— (l—e“)e S e-i20e: ny even. (7) 
a i=0 
The summations in equations (6) and (7) are the sum of (n+1) /2 
and n/2 terms respectively of the geometric progression of common 
ratio e***. Equations (6) and (7) thus become 


(m+1) 


AE 1—e 
e(nc) =— [ht en ret n odd (8) 
a 


and 


AE tee 
2(7ne) =——— (1b — ¢-*) —______¢€“; neven. (9) 
a sip eee oes 

When v is odd e(ne), the value of « at the end of 7 half periods, is 
also the value of « at the end of (n+1)/2 “stimulatory periods” 
[intervals during which E(t) = EF]. When 7 is even «(nc) is the 
value of « at the end of n/2 “rest periods” [intervals during which 
E(t) = 0]. Thus (8) or (9) gives « directly in terms of the number 
of stimulatory periods or rest periods when the total time interval 
terminates with a stimulatory or rest period. It is also clear that if 
m is an odd integer then, from (8) and (9), 


e([m + 1]c) =e(me)e, (10) 


the physical interpretation of which is quite apparent. 
When F(t) has the form 


BGy— fh, 7h) Stan lit nk 


=O) ea teantn = t <1) (A) (11) 
n=0,1,2,---: 
equation (3) gives for « at the end of the nth ¢,-period 
n-1 (k+1) tatkte 
e(nt, + [n—1]t.) =AED Pela Cae Lg Hl (12) 


k=0 K(titte) 
When the definite integrals in (12) are evaluated the result is 


AE ir 
e(nt,'+ [n—1]é.) =— (1—e*) De" 
a i=0 
(13) 
A E 1 a ener 
=— (1—e*) 
a 1 — e-4t 
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where + = t,+t.. In quite the same manner the value of « at the 
end of the nth ¢.-period is found to be 


AE 1—e*4 
e(nt, + nt.)=— (1 — e*%) es ———. (14) 
a i tse 
Equations (13) and (14) are identical with the results of Rashevsky 
(1948). They yield the physically obvious relation, equivalent to 
(10), 


e(nt, + nt.) =e(nt, + [n—1]h)e™. (15) 


When t, = t2 = c, equations (8) and (9) are recovered from equa- 
tions (18) and (14) respectively. 

By the same general method two further results which are of 
possible interest may be readily obtained. Let H(t) be given by 


E(t) =E sin wt, n(& +/o) <tsnt,+ (n+1)2n/@ 
=0, Nt, + (n+ 1)a/o <t< (n +1) (4.4 2/m)_ (16) 
R= Openers 
the value of « at the end of s stimulatory periods is found to be 
&(Sr wrt. [Ss 11) 
AE ow ** 


Sao oe Dd eB(r/o+ta) [qa sin a — w cos al (17) 
a OO” k=0 


X [cos(k + 1)a— e*/* cos ka], 


where 8=s—k—1 anda= kot.~ 
If t2 = a/w then (16) gives the half wave rectification of E sin 
wt and (17) becomes 


A @ %-1 
6 (128 — 1] nfo) oe aria tee 

a? ae a” i=0 

AEo 


1 
= a! = E-27/w) ; 
a? + @? 1 — e2t4/o 


(18) 


— e-8:2aT/w 


Since S stimulatory periods correspond to a total of n = 2s — 1 half 
periods, equation (18) with 5s = (n+1)/2 gives e(nxa/w), the value 
of « at the end of an odd number, 7, of half periods. 
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ON PERMEABILITY DETERMINATION OF 
THE ERYTHROCYTE 
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COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


A procedure is indicated to test the assumptions underlying the 
theory of determination of the permeability coefficients of the erythro- 
cyte from its volume changes. 


A method for the determination of the permeability coefficients 
of the erythrocyte from its volume changes caused by the penetra- 
tion of an isosmotic solution into the cell is well known. The differ- 
ential equations which are encountered in that method are non-lin- 
ear. M. H. Jacobs (1934), who developed the method, gave numeri- 
cal tables for the solution of these equations. It is the purpose of 
this note to indicate a procedure for checking the assumptions in- 
volved in the theory of this method, for instance, the assumption of 
independence of the permeability coefficients on the solute concen- 
tration. E Schiddt (1933) has shown that a sufficient number of 
measurements can be carried out to represent the relation between 
the erythrocyte volume V and the duration of diffusion t by a smooth 
curve V(t). The availability of such a smoothed relation is required 
for the application of the procedure here indicated. 

Let us consider an erythrocyte surrounded by an isosmotic solu- 
tion of a particular solvent. Let S; be the number of solute mole- 
cules inside the cell at any time ¢ and W, the number of solvent mole- 
cules therein. Let [S;] [Sz] [W:] [Wz] be the concentrations of 
those substances inside and outside of the cell respectively expressed 
as number of molecules per cm’. Let [o,] be the concentration of 
other osmotically active substances inside the cell whose diffusion 
through the cell membrane is assumed to be negligible, so that their 
amount o; inside is constant. Let hg and hw be the permeability co- 
efficients of the solute and of the solvent, and A the area of the 
erythrocyte. The definition of the permeability coefficient gives the 
following differential equations: 
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dS,/dt = hgA ({Sx] — [S1]), (1) 
dW,/dt = hwA({Wa] — [W7] Js (2) 


The particular shape of the erythrocyte makes it possible to assume 
that within the range of swelling considered the cell maintains its 
area constant (cf. Ponder, 1948, p. 82). This implies that the vol- 
ume change of the cell is unaccompanied by changes in the stress of 
the cell membrane, which requires a constancy of the external and 
internal pressures during the process of diffusion. Consequently, 
assuming that the relationship between the pressure and the con- 
centration is ruled by the law of perfect gases, we have: 


[S,] + (W.] + [or] = [Se] + [We] =C=constant . (3) 

The experimental conditions are, in general, such that [Sz] and 
[W,] can be considered individually also as constant. Then: 

dV /dt = Cd(S;+ W, + 0;)/dt = (dW,/dt)/C, (4) 


if we neglect the change of the amount of the solute with respect to 
that of the solvent. Substituting equation (4) in equation (2) and 
eliminating [W;] and [W;] by means of equation (3) we obtain 


ds/dt = [1— (s/v)]o, (5) 
dv/dt=[(s + 1)/v—1]o, (6) 
where 
$s = S,/or oO = [Sz] V/or oo hsA [Sz] /ox gb == hwA [Sz]?/ (oC). 


Equations (5) and (6) are Jacobs’-equations. An initial condition is 
(0 j= 0 » Le., no solute exists inside the cell at the beginning. With 
these conditions equations (5) and (6) can be written in the form: 


t 


v dv 
s(t) sven [ edt —=——+y—1, (7) 
x dt 


@ 


t 
where I(t) = { dt/v(t). Since w is known from measurements on 
0 


penetration of the pure solvent into the cell (cf. J acobs, loc cit.) and 
v(t) is known from the measured changes of the erythrocyte volume 
in the presence of the solution, all quantities in the middle and right- 
hand sides of equation (7 ) are numerically known, except «. Thus 
the problem of finding o is reduced to that of fitting a given curve 
R(t) to a one-parametric equation 
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o[J (t)]-? | [J (t) ]° dt = R(t) (8) 


0 
where FR and J are known functions, R being the right-hand side of 
equation (7) and J(t) = exp [J(t)]. A value of o is obtainable from 
Jacobs’ tables, and the difference between the two sides of equation 
(8) gives an indication as to the accuracy of the fit and the correct- 
ness of the assumptions underlying the theory. Equation (8) can be 
also used to determine oc. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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ON THE DIFFUSION OF METABOLIC INTERMEDIATES 
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A simple model of diffusible metabolic intermediates, produced in a 
spherical cell, is postulated. The rate of production of the intermediate 
is taken to be a linear function of its internal concentration. The inter- 
mediate is produced by the degradation of some precursor substance 
at a rate which is assumed to decrezse with increasing internal concen- 
tration. Assuming the medium in which the cell is enclosed to be infinite 
in extent, and the external concentration of the intermediate constant, 
the time dependent and steady-state solutions of the diffusion problem 
involved are obtained. It is shown that production of the intermediate 
will not cease, for any value of r < r,, where 1, is the radius of the 


cell. Possible application of this approach to the study of the occur- 
rence of the lag phase in bacterial growth processes is indicated. 


The elucidation of the réle of diffusible metabolic intermediates 
in the overall activities of the metabolizing cell is of fundamental im- 
port in cellular biology. The studies of C. N. Hinshelwood (1946) on 
the kinetic aspects of bacterial growth have been concerned, in part, 
with the occurrence of the lag phase and the diffusion of certain 
essential intermediates. He views the duration of the lag phase as 
depending on the presence of a certain amount of the intermediate 
in the culture medium when a transfer is made from a growing cul- 
ture to a new medium. 

If we consider the phenomenon of lag as being due to a neces- 
sary readjustment in the metabolic processes of the organism that 
must occur before the organism can enter the logarithmic growth 
phase, then there may exist a relationship between the rate of a con- 
secutive set of reactions taking place within the cell and the rate of 
change of the external and internal concentrations of some interme- 
diate produced. 

The mathematical biophysics approach to the problem of diffu- 
sion of metabolic intermediates will be to postulate certain metabolic 
models that exhibit different modes of intermediate production, and 
to examine these models for the relationships between internal and 
external concentrations of the intermediates, and the influence they 
have on the behavior of the metabolic model. 
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DIFFUSION OF METABOLIC INTERMEDIATES 
Sa aae 
---»Ss, —>s, ——> + 


So 
FIGURE 1 


In this paper we will consider a simple model such as depicted 
in Figure 1. Consider a spherical metabolizing cell of radius 7) in 
a medium of infinite extent. Let a simple reaction, such as shown 
above, be taking place where a precursor substance S, gives rise to 
another substance S.. The substance S, is degraded in the next step 
of the reaction to give a substance S; and an intermediate which we 
will designate S,. We will denote by c; the internal concentration 
of the intermediate, and by c; its external concentration. 

Consider the general diffusion equation (Rashevsky, 1948) 


0G; 
ot 


where D; is the internal diffusion coefficient, c; the internal concen- 


tration, and q the rate of production of the intermediate. Let us 
assume that 


D; V?¢q+q= ’ (1) 


q=a—be;,, (2) 


i.e., the rate of production is a linear function of the concentration, 
varying from a when none is present to zero when c; = a/b. 

The symmetry and homogeneity of the spherical cell indicate 
that the concentration c; will be a function of ¢ and the single space 
variable r , the distance from the center of the cell. In spherical co- 
ordinates equations (1) and (2) reduce to 


Ores Pay BY (Of 0G; 
D; + — + (a—be;) = sare (3) 


0 r? Tr Or 


This partial differential equation holds within the cell. 
The rate of flow ®, of the intermediate per unit area across the 
boundary of the cell, at r = 7, is assumed to be proportional to 


A4nr,?(¢; — c;), thus 
©, = Anr,*h (¢; ssa C;) ; 


where h is the permeability. Across any spherical shell of thickness 
Ar wihin the cell we expect the intermediate to diffuse at the rate 
Anr,?(Ac;/Ar), or across any spherical surface, 
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6) Ci 
®, =— 4anr,?D; : 
One j. 


But since the intermediate must be conserved, we have at the bound- 
ary 


&, = , , 


and so 


/ Os 
ee ett (4) 
OF: 


Let us donate by ¢, the initial concentration of the intermediate with- 
in the cell and assume a constant outside concentration. In con- 
sidering the external concentration to be constant we assume the 
medium to be well-stirred or the external diffusion coefficient D; = o. 


Solution of Problem. If we make the substitutions 
v = rD;"(a— be;) 
and 7 = D;t , the boundary value problem becomes 


ov b ov 
Snel va (5) 


with the conditions 


Ov 
=) = Ghrhe (ris) (6) 
or iy 
and 
vi 
v(r, 0) =— (a— ben). (7) 
D; 
The additional condition : 

v(0,7) =0 (8) 
arises from our definition of v(7, 7). In (6) we have used the con- 
stants 

Toh ; Wa D,— hry 
Sener (a— Cy) ; <9 Da ° 
Now let 


v(r,7) =wtr,r) +27), (9) 
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where 
2 Ow 
UE 3a ee (10) 
or? D; Or 
and 
d?z b 
—_— —z=0. (Eby 
dr? 2D; 
Further, require that z(0) = 0, and 
dz 
— = Go iZteye 
dr 


while w(7, 7) satisfies the conditions 

ow 

ee Hite) (12) 
and 

w(0,7) =0. (13) 
We see that 
G sinh rk 

[k cosh 7. —H sinh 7k] 


satisfies the conditions imposed on z, and we now find from (7) and 
(9) that 


2(r) = 


7 G sinh rk 
wir, 0) aiece (a — be,) — (14) 


i [k cosh ».k — H sinh ron] ; 
A b 
In the above expressions we made the substitution k = J ae Equa- 


tions (10), (12), (13), and (14) constitute a boundary value prob- 
lem in w(r, +) which is to be solved. 
Let 


w(r,7) =R(r)T (zr), 


and substitute in (10), writing it in the form 


=— + he. (15) 
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Since the function on the left can vary only with r and the one on 
the right only with 7 , they must both equal a constant, say —A?. From 
(15) we have the two ordinary differential equations 


Rk’ +vVR=0, (16) 
and 
Deemer A), (17) 
From (12) we find that 
(R'— HRY, =0, (18) 
and from (13) 
R(0) =0. (19) 
The system (16), (18), and (19) has solutions 
R(r) =E, sin apr, (20) 
where j,(n = 1,2, 8.----) are roots of the characteristic equation 
Aw = 11 - tan Ago « (21) 


Roots of this equation can be determined graphically, or found in 
mathematical tables (Jahnke and Emde, 1945). 
Equation (17) has the solutions 
T, = Fye-Fn)7, (22) 
therefore the product 
hele =i Koes sin or (23) 


satisfies all requirements for w(r, 7) except (14). 
We now write 
a0 Ti ea nl SiN Aals 
n=1 
and since the functions sin 4,7 are orthogonal in the interval0 <r < 7% 
(Sommerfeld, 1949) we choose K, as the Fourier coefficient. in the 
series from w(r, 0), as given by (14). Thus 


wir, 0) = > KK, Sin Anis 


n=1 
and 
4Ay es 


K,= : 
{2An?o — SiN ZAnTo} © 


w(r,0) sin dA,rdr; 


hence 
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Ah, (a a= bCo) 

| eT ae 
[2An%o — Sin ZAnto] | Di An? 

G(k sin Ano cosh kro — An Sinh kro cos AnYo) 

ee Soe cee ee ee ee 


(k2 + An?) (k cosh kro — H sinh kro) 


(sin AnTo — AnTo COS AnTo) 


The solution, ¢;(7, t), can now be written 


a Di * D/Dit\n2)t ai 
e,(r, t) =——— > Kye 0)? sin An? 
b a n=1 


(24) 
Gd; sinh kr 
x a 


+ b{k cosh kr, — H sinh kro} T 


The Steady-state. From the time dependent solution we find 
that in the steady-state 


a hr,? (a — be;) sinh kr 
SS ~ 2(25) 
b b{Dirok cosh kr, — (Di — hro) sinh kro} r 


We know from our definition of qg as a function of ¢; that if 
c; = a/b the intermediate will no longer be produced. Hence we are 
interested in determining if for any r < 7, the internal concentra- 
tion will equal the above value. If the second term on the right in 
(25) remains positive c; will, for all values of r < 7, be less than 
a/b. Consider the denominator of the coefficient of (sinh kr)/r: 


P(%,) = b{Dirok:cosh kro — (D; — hry) sinh kro} . 
Let us rewrite P(7,) as follows: 
P (70) = b{D; (kro cosh kr, — sinh kro) + hry sinh kro} . 


Both sinh kr, and cosh kr, are monotonically increasing functions. 
The expression kr, cosh kr, — sinh kr, is zero for 7% = 0, and its 


derivative k?r, sinh kr, is positive everywhere. Hence P(7%) is always. 
positive. 


At the center of the cell, 7 = 0 , we have 


: a 
o—lime,=— 


r>0 


khr,? (a — be;) sat 


b{D; (kr, cosh kr, — sinh kr.) + hry sinh kro} ; 
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The distribution of the internal and external concentrations is the 
system shown in Figure 2. 


Yp 


: | 


) 


| 


° fy 7 —>» 


FIGURE 2 


The application of this approach to chains of autocatalytic re- 
actions and cyclic reactions is now being considered, with the hopes 
that it may be useful in the study of the role of diffusible metabolic 
intermediates in growth processes. 

The author is indebted to Dr. H. G. Landau for reading the 
manuscript and for several helpful suggestions during the course of 
preparation. 

This work was aided in part bya grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of The University of Chicago. 
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A NOTE ON RAPOPORT’S APPROXIMATE FORMULA FOR 
THE INPUT-OUTPUT CURVE IN THE CASE OF 
BILATERAL PRE-INHIBITION ‘ 


S. LEVINSOHN 
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The accuracy of an approximation method used in a certain input— 
output problem of randomized stimuli is evaluated. The curves derived 
trom it are shown to be close approximations to those derived by a sta- 
tistically ‘exact’? method. 


A. Rapoport (1950) derives the following expression for the 
output of a neuron stimulated by two other neurons in the case of 
“bilateral pre-inhibition.” 

20 (eo 1) 
<< (1) 
2x? 6? + 4% 6 + €% 


where x is the frequency of the outside (Poisson-randomized) input 


and x; the output of the neuron. 
By an approximate method, Rapoport derives the expression 
2% 
=, (2) 
2x O:taL 


To determine the maximum error obtained by working with 7; 
instead of with x, we can proceed as follows. Let #;/x = R. Then 
nae aa tote 


R=1+————_-. (3) 
Qu? 0? + 3204+ 1 


Let x 6 = y; then 


dy 
(2y? + 8y + 1) (1—e”) — (yt+e"—1) (4y + 38) 
Bete Soe 
39 


(4) 
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To find the values of y for which R is a maximum, we set R’=0. 
Then 


—2y? + dy + 4—e4(2y? i+ Ty + 4) = 0; (5) 
—2y?+4y+4 
Speen (6) 
2y? + Ty + 4 
Let 
—2y? + 4y +4 
——_—_———_——- =r. 
2y? shal yined 


Aty=0,r=1,andR=1 (aminimum value of Ff). 
Likewise for y= », R= 1. Calculation of r and e~ to three 
decimal places gives us the following table: 


Yy f ev 
0.00 1.000 1.000 
0.01 0.992 0.990 
0.10 0.928 0.905 
0.50 0.687 0.606 
1.00 0.461 0.368 
2.00 0.154 0.185 
2.10 0.130 0.122 
2.15 0.118 0.116 
2.16 0.116 0.115 
Pest 0.114 0.114 
2.20 0.107 0.111 


‘The maximum error in the use of the approximate formula oc- 
curs when r = e”, ie, at y= 2.17. Then R = #3/x; = 1.0759. Thus 
the maximum error is seen not to exceed 8%. 
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ON THE MATHEMATICAL THEORIES OF EXCITATION 
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The question of constancy of the velocity of propagation of exci- 
tation is discussed on the basis of the two-factor theory, and differences 
between this and Blair’s theory are examined. 


The mathematical theories of excitation express the effect of a 
stimulus as the magnitude of a quantity, ¢ in Blair’s theory, with 


de/dt=cal—pe, (1) 
and « — 7 in Rashevsky’s theory, with 
a e/dt =oae 1 — pe(e— &) 


dj fdt oid — pi(7— Jo). 
The o’s and p’s are constants characterizing the biological system, 


t is the time, J the magnitude of the stimulus, the initial conditions 
being <«(0) = 0 for equation (1) and 


(0) =e ,9(0) — 4. with Jo > & (3) 


for the system (2). In Blair’s theory excitation occurs at a moment 
r at which ¢e reaches a certain magnitude h which is also considered 
as a characteristic of the biological system. Thus, from equation 
(1), the response time 7 is determined by the relation 


(2) 


ot er { ser Cur= 1 where) ot —=c/h. (4) 
0 


In Rashevsky’s theory the excitation occurs when « = j, and the 
relation determining the response time 7 can be put in a form simi- 


lar to equation (4). In fact, putting «= « — =, 7 = Jj — jo, and 
omitting the bars for simplicity, equations (2)—(3) become 


Al 


42 MATHEMATICAL THEORIES OF EXCITATION 


de/dt=ocl—peé, (0) =0, 
| coe (5) 
| dj/dt=o;l—pii, (0) =0, 


and the condition of excitation « = 7 becomes « — 7 = to Eo =H. 
Thus, through the use of equation (5), the response time 7z is deter- 
mined in Rashevsky’s theory by the relation: 


Ge eX (= pa) [exp (ret) Oat 
(6) 
a; exp(—p7) | exp (p;t)I(t)dt=1, 


where oc* = oc/H , o;* = 0;/H. It is known that equation (6) is also 
obtainable from Hill’s theory. Blair’s theory involves two independ- 
ent parameters (o*, p) and Rashevsky’s theory involves four (ce’, pe , 
c;", pj). Both theories consist of associating a time 7 with each func- 


tion I(t); equation (4) or (6) determines 7 as a functional of the 
stimulus I(t). 


The propagation of excitation is examined in this paper on the 
basis of equation (6) in a similar manner as it has been done on the 
basis of Blair’s theory (Opatowski, 1950). Following Rashevsky’s 
treatment we consider a nerve fiber to which a constant current I, 
is applied at one end, x = 0. Assuming the fiber to be of infinite 
length, the longitudinal current at a distance x from x = 0 is, in the 
steady state, 


LJ eo", (7) 


where f is a constant dependent on geometric and electric charac- 
teristics of the fiber. If a possible contribution of internal energy 
sources to the process of excitation is neglected, the substitution of 


equation (7) into (6) determines a function z(r) which describes 
the propagation of excitation: 


(oe*/pe) [1 — exp(— pet) | 


8 
— (0;*/p;) [1 — exp (—p;r) ] = e8*/I,. © 


The velocity of propagation of excitation dx/dz is obtained by dif- 
ferentiating equation (8): 


di:/dx = Iye-*[o¢" exp (— per) — oj" exp(— pyt) 1/8. (9) 


ioe is eliminated here by means of equation (8) an explicit ex- 
pression of the time dependence of this velocity is obtained: 
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, (10) 
ge €XP(— pet) — o;* exp (— pit) 


Me ye lea oexd Cae ia explean 


A comparison of equations (6) and (4) shows that if pe = pj 
the two-factor theory reduces to the one-factor theory with p= pe=p} 
and o* = o.* — o;*. Therefore if the two-factor theory is to yield 
more than a mere correction of the results of the one-factor theory, 
the constants pe and p; must not only be different but also of differ- 
ent orders of magnitude. As a matter of fact, this is the actual sit- 
uation in the experiments analyzed by A. V. Hill et al in which p; 
is, in general, 10 to 200 times smaller than p-; in one case only (frog 
sciatic nerve) this ratio reaches a value as low as 1.6. 

Since there are experimental indications that the velocity of ex- 
citation does not vary greatly during its propagation it is of interest 
to discuss such approximate independence of dx/dr from 7, on the 
basis of equation (10). Since the right-hand side of this equation 
cannot be made exactly independent of 7 for any choice of the con- 
stants, the velocity of propagation cannot be exactly constant in a 
mathematical sense. 

The velocity of propagation of excitation at x = 0 is obtained 
from equation (9) by putting = 0 and 7 = 7, where 7, is the re- 
sponse time at x = 0. A necessary condition for the excitation to 
propagate at all is dx/dr > 0 atx =0. Therefore, if the right-hand 
side of equation (10) has a positive zero, r = 7, the latter must be 
larger than a possible duration of propagation of excitation along a 
whole nerve because otherwise dx/dr could not be approximately con- 
stant. But from the numerator of the right-hand side of equation 
(10) we have: 


— (p; ore pe) tn (a;*/oe*) F 
Thus even if p; << pe, so that p; — pe is numerically large, a choice of 
a sufficiently small o;* with respect to ce* will make 7 as large as 
desired. That such a choice of constants in principle does not con- 
tradict equation (8) can be seen in the following manner. Since, 
according to our discussion, p;/pe = 7 and o;*/o," = s have to be as- 
sumed small numbers, we have from equation (8): 


pe 0°*/(oe* I.) = 1 — exp(— per) 
— (s/r) [1— exp(—7 per) ] ~ 1— exp(— pet)— 5S pet, 
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which is positive because of the smallness of s, as it should be. Thus 
an approximate constancy of the velocity of propagation of excita- 
tion seems to be compatible with the two-factor theory even af one 
neglects any possible contribution of internal energy sources to the 
process of excitation. Lack of sufficient information on the magni- 
tude of the constants appearing in the theory makes it impossible to 
arrive at more definite conclusions. 

N. Rashevsky assumed that the mathematical expression of the 
distribution of the current during excitation was obtained by a rigid 
displacement of function (7) along the «z-axis with « = 0 always at 
that point of the fiber which was excited at the moment considered. 
This can be looked upon as a way of taking into account the contri- 
bution of internal energies to the process of propagation of excita- 
tion. It can be easily seen that such a theory cannot lead to a mathe- 
matically constant velocity of conduction of excitation, although it 
is compatible with an approximately constant velocity as has been 
shown by Rashevsky. 

In fact, if E(t) is the distance covered by the excitation during 
time t, equation (7) must be replaced now by: 


T= I, eb bs@, (11) 
Putting, for 7 =e, 7: 


t 
Jn =I on EXP (— pnt — 6 x) { exp [pal +e? a dh 


we obtain from equation (5): 
saedi= Je lth (12) 


By an argument similar to that used by the author in connection 
with Blair’s theory (Opatowski, 1950), the velocity of propaga- 
tion of excitation is constant if, and only if, [0(e — 7) /Ox]ex. and 
Lo(e — 7) /ot]exe , calculated from equation (12), are constant, where 
the subscript exc indicates that the functions in the brackets are 
to be computed at that time t at which the point x becomes excited, 
the condition of excitation being e — 7 = H. Since [elec seu) 
it is seen from equation (12) that 


Lee 1) /02] exe = — 8 H = constant, 
[a(e— 9) fet] ecole las michy) se Lp; Jj — pe J aidecé . 


But [Je — JyJexe = le — jl exe = H; consequently the right-hand side 
of the last equation is constant if, and only if, p; = pe, but this can- 
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not be accepted because it would reduce the two-factor theory to the 
one-factor theory. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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A NOTE ON FACILITATION AND THRESHOLD PHENOMENA 


NORMAN ANDERSON AND ANATOL RAPOPORT 
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A neuron subjected to a Poisson shower of stimuli responds only 
if h stimuli impinge upon it within the time interval o. It is shown 
that the derivative of the input-output curve cannot exceed unity. 


One of the authors (Rapoport, 1950) derived the input-output 
curve for a single neuron under the following conditions: 


1. The neuron is subject to a Poisson shower of stimuli. 

2. The neuron responds only to the hth stimulus of a group 
which falls within the interval o , provided none of the h stimuli im- 
pinges during the refractory period 6. 


The expression for the input-output curve turned out to be 


f(%,o,6,h) 
(1) 
x[1—e°? B,2(% o)] 


< x[1—e*? E,.(4@0o)) +1+ (h—1) [1—e*? Ena (2 aye 


k 
where E;,(z) = > 2//7!. 

j=0 
The shape of the curve is sigmoid (cf. Rapoport, 1950, Fig. 3). A 
question naturally arises whether by a proper choice of the para- 
meters, ¢, 6, and h, the curve can be made to exhibit a “quasi- 
threshold” effect, i.e., a sudden rise from a very small value to prac- 
tically its asymptotic value. Such a “threshold” effect is certainly 
present if the stimuli are regularly spaced, as has been shown in the 
above-mentioned paper. It was conjectured, however, that for the 
case of a Poisson shower (randomized stimuli) no such threshhold 
effect could be obtained by any choice of parameters. This is in con- 
sequence of the fact that the derivative of/ox of the input-output 


AT 
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curve (1) is bounded for all values of s,o, 0, and h. This property 


will now be proved. 
In expression (1) call the threshold h, instead of h and let 


h=h,—2. We also introduce the following notations 


: 
2 (ga)? 
on (2 0) =1—e* Ey (ao) =e D , (2) 
h+1 t! 
D=éa[l—er kh, (#0) ] + Pa WA) ee 
Zz (3) 
——O 2 Gyr lat (i eke 
With this notation, 
Xv Ap 
(4) 


ee 
OL Onan ik ae (h + 1) Gis 


Let primed quantities represent derivatives with respect to x. Then 


1 
f= 


De (ay, Se 45 a,)D—x An D | 


1 
=5 (a, + % az’) (3°xa; 1S (h + 1) ojs1) — 2 a; (6 a, + 0 X ay’ 
ae + Det) | (5) 
1 
ant 6 % an? + a, + (h +1) Gi Gh + 0 UG? a, a, + £0, 


+ (h +1) a og 8 toe = 1) ee «| 


Now 
2 (ua)? ~ (xo) k 
Oy i es 0 62° bey) 2 8 eh ’ (6) 
k+1 th k t! k! 
and 
( k+1 


Dy = (k +1) ee = (k 1) (Oe+1 — Oxs2) . (7) 


(k+1)! 


Thus we see that 
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1 
r=z [nc Sia (h =a 1) ops) — (h Se 1) (h AF 2) on (Onr2 — Ansa) 
(8) 
Tait er) pen —— One C4 (Pe 1.) aja) | ‘ 


We have 0 < an < a, < 1 forall &k. Also 
i= 6 40) Fd (h+1)aan 214+ (h + 1) ann. (9) 


Hence, 


ics : 
s= 


tS (hee ews |] + (h +.1) (aia — ama) | 


_ Lt (ht Vana] Len + (e+ 1) (ann = an) ] 
[1+ (h + 1) ain]? 

__ Oo + (h +1) (ena — ame) 

Mee tet (hisst) oe | 


(10) 


Then 7’ < 1, since 
On —- (h =r 1) (Ona — Ons) <1 
1+ (h +1) Gin 


This last inequality is true because a, < 1 and (h + 1) (on — ane) 
< (h a 1) Oss ° 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of The University of Chicago. 
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SUGGESTIONS FOR A MATHEMATICAL BIOLOGY 
OF SOME CULTURAL DEVELOPMENTS 


N. RASHEVSKY 
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In connection with the author’s previous studies on the effects of 
imitation in social behavior it is shown how, owing to such effects, two 
societies which are both characterized by the same distribution functions 
for different abilities and tastes may differ very greatly in their respec- 
tive outputs of scientific and inventive work. The course of development 
which occurs in a given society may be determined by purely accidental 
initial conditions. A theory of quantitative relations between some cul- 
tural and socioeconomic quantities is suggested. ae 


It is now generally accepted that human races do not differ ap- 
preciably from each other with respect to their mental abilities. In- 
asmuch as within any given race different mental abilities vary from 
individual to individual and are distributed within the population 
according to some distribution function a precise comparison of dif- 
ferent races requires a comparison of those distribution curves. 
While such comparisons are, perhaps, not available with sufficient 
precision, it is plausible to assume that those distribution curves are 
the same for different races. But if so, the incidence of genius in a 
given field should be the same for all races. ‘The circumstance that 
this is patently not the case must then be explained as being due to 
cultural, learned conditions, as well as to differences in environment. 
This is plausible enough, and may be considered as a good qualitative 
theory. As always, however, a real check by observation is possible 
only with a quantitative theory. In this paper; we shall outline a 
quantitative approach to the problem. If developed, this approach 
should enable us not only to understand the reasons for differences 
in certain types of achievement in different races, but should also 
lead to quantitative relations which are verifiable by,,observation. 

In chapter xii of our book (Rashevsky, 1950). we. have consid- 
ered the effects of imitation upon the behavior of a social group in 
which each individual can exhibit one of two mutually exclusive be- 
haviors, R, or R.. The tendency of each individual toward either 
R, or R. is measured by a neurobiophysical quantity ¢. Positive 
values of 6 mean a preference for R;. Negative values mean a pref- 
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erence for R,. The larger the absolute value of ¢, the greater the 
corresponding preference. An individual with ¢ = 0 has no pref- 
erence for either FR, or FR, . 

We have seen that, even if the distribution function N(¢) of the 
quantity ¢ in the population is symmetric with respect to ¢ ie 0 so 
that on the average there is no preference for either ie Oreies in the 
population as a whole, nevertheless, owing to imitation, the majority 
of the population will exhibit either RF, or Fz , depending upon purely 
accidental initial conditions. 

The behaviors R, and R, may be of any kind. Thus R, may stand 
for a particular fashion in dress, R, for a different fashion. Or_ Fk, 
may stand for militaristic, R. for pacifistic behavior. We may, there- 
fore, also consider R, as representing a behavior which consists of 
the avoidance of any except the most necessary activity. The Italians 
call it dolce far niente. In common English parlance it is known as 
laziness. The opposite behavior to that, namely R., is an activity 
sometimes beyond the bounds required for the satisfaction of our 
daily needs. Such men as scientists, inventors, and artists exhibit 
for the most part behavior R,. It is legitimate to apply the conclu- 
sions of loc. cit. to the case of those two mutually exclusive behav- 
iors. There hardly can be any doubt that the attitude of an indi- 
vidual toward work is determined not only by his innate laziness or 
industriousness, in other words his ¢, but also by the behavior of 
the other members of the society in which he lives (the quantity y 
of loc. cit.). 

Now let us consider the following theoretical case. Let us take 
two societies, I and II, which have the same distribution function 
N(¢). Let us further consider, for simplicity, the case in which 
N(¢) is symmetric with respect to ¢ =0. However, let society I be 
in a configuration in which R, is preponderant, while in society II R, 
is preponderant. We shall speak of configurations R, and R. corre- 
spondingly. As in loc. cit. we shall denote by X, the number of indi- 
viduals which exhibits R, in society I, in the stable configuration, 
while Y, shall denote the number of individuals which exhibits R, in 
society I. For society II we shall denote the corresponding quantities 
by X. and Y,. Let N, and N, be the number of individuals in the two 


cei Because of the symmetry of N (¢), it follows from loc. cit 
a 


peered (1) 
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ser epee oe ee (2) 


The X, individuals in society I will do only the minimum amount 
necessary to keep them alive, while the Y, individuals will apply 
themselves to all kinds of different activities above and beyond their 
minimum needs. Under rather general conditions, as we saw in loc. 
Cit. 

L—a<<il. (3) 


Thus society I will consist mainly of poor individuals without any 
interests outside of their minimum daily requirements. On the other 
hand, society II will produce considerably more than its necessary 
average minimum subsistence. 


As was said above, men like scientists and inventors exhibit be- 
havior FR, , and hence have a large negative ¢. But a large negative 
¢ is certainly not enough to make an individual an inventor or scien- 
tist. He also needs sufficiently high ability. While scientific and in- 
ventive ability are probably of a rather complex nature and involve 
a large number of factors, we shall, for simplicity, consider here the 
theoretical case of only one ability a, distributed according to a dis- 
tribution function V (a), such that 


{ Vi@da=1. (4) 


Let both societies have the same V (a). 

But a large a does not insure that a man will become a scientist, 
inventor, or artist. He must have both a large a and a large nega- 
tive ¢ + w, that is, he must have a large a and at the same time ex- 
hibit behavior R.. There are individuals with undoubted talent who 
never achieve anything because of behavior R,. The classical Rus- 
sian literature is full of such types, the most typical being Goncha- 
roft’s Oblomoff. 


If we arbitrarily consider an individual as outstanding or a 
genius if his a exceeds a prescribed value a*, then, putting 


[/ viaaa=s, (5) 


we find from (2) and (5) that the percentage G, of scientists and 
inventors in society I is equal to 


G,= (1—a)6, (6) 
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while the percentage G, of scientists and inventors in society II is 
given by 


G,=afp. (7) 
In the case in which inequality (3) holds, we have 
Gi<<G. (8) 
The value of 1 — a may be as small as .01 or even less. In that case 
a > 10°. 


‘1 
The actual situation is much more complex. We shall now con- 
sider a slightly more complex and realistic case. 


As we have seen, society II produces more goods and commod- 
ities of all kinds than society I. In monetary terms we may say that 
its average per capita income i is larger. The income of society I 
will be distributed according to some distribution function U,(%), 
such that 


i U,(i)di=1. (9) 


The income of society II will be distributed according to a corre- 
sponding function U.(i). Several possibilities may be considered. 
All individuals who exhibit R, may belong to the higher income 
bracket. This is a rather natural assumption. Or this may not be 
strictly the case, but negative values of ¢ + w and income may be 
positively correlated. Finally, we may consider the case in which 
there is no correlation between ¢ and 7. We shall consider here only 
the last case, without prejudice to the other two. The case is not so 
unnatural as it may first seem. While, say, in society I Y, individ- 
uals exhibit R, at any given time, the group Y, is not composed of 
the same individuals at different times. Moreover, as we saw in chap- 
ter ix of loc. cit., pure chance may perhaps play a role in the distribu- 
tion of wealth. 


Scientific, artistic, and inventive activities are, as a rule, not 
remunerative per se. While in some contemporary societies some in- 
ventors have become wealthy through their inventions, the opposite 
seems to have been the rule in the course of even recent history. Max- 
well probably never dreamed of the wealth which his discoveries 


would bring half a century later to General Electric, Westinghouse, 
and the Edison Companies. 
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A non-remunerative occupation can be pursued essentially only 
in hours of leisure, when one is free from the toil for survival. With- 
in limits we consider that the greater the income of an individual, 
the more leisure time he has. Perhaps it is better to consider in this 
connection accumulated wealth w rather than income. A business- 
man who has a large income primarily because he bends all his ef- 
forts to making money may have no leisure time at all. A man who 
inherits wealth may have plenty of leisure time. An individual who 
has to work 12 hours a day to make a bare living does not have time 
to meditate on the Theory of Relativity. Thus, in general, the greater 
w, the greater the available fraction 6 of time which an individual 
may devote to science, art, or invention. Inasmuch as 6 < 1, the 6, 
w curve must approach 1 asymptotically as w increases. 

The results of any scientific or inventive activity depend on the 
one hand on the ability of the individual, on the other on the amount 
of time he puts into this activity. The greater 6, the less may bea, 
in order to achieve definite valuable results. Since we consider here 
the case in which 6 is a monotonic increasing function of w, we may 
say that in a group of individuals with a given w only those will con- 
tribute to science and invention for whom 


a>u(w), (10) 


where u(w) is a monotone decreasing function of w. It may, for 
example, be of the form 


i 
u(w)=—+a'; (11) 
w 


where a* has the same meaning as in equation (5). Expression (11), 
together with (10), implies that even for a very high w, when @ is 
practically equal to 1, a minimum ability a* is necessary to make a 
good scientist or inventor. 


Consider now society I. With the assumptions made above the 
distribution of wealth among the Y, individuals is given by the same 
distribution function as the distribution of wealth in the whole so- 
ciety. Actually the distribution curves of income and wealth are very 
similar and we may, therefore, assume that w is distributed in the 
same way as isi. In a group with a given w the fraction of scien- 
tists or inventors is given by 


° V(a)da=T(w). (12) 


u(w) 
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Hence the fraction of individuals who have wealth between w and 
w + dw and who are scientists or inventors is equal to 


(1—a)T(w) U,(w) dw, (13) 
and the total fraction in the whole population at all levels of wealth 
is 

G,= (1—a) f T (w) Us (w) dw. (14) 
0 
Similarly, for society II we find 
G.=a i T (w) U2 (w) dw (15) 
0 

Denoting by T, and T, the mean values of 7'(w) in the two societies, 

we may write equations (14) and (15) thus 
G,= (1—a)T); (16) 
Ge —a fe . (17) 


If, as would naturally be the case, V(a) is zero for a = 0, has a 
maximum for a positive value a, of a, and tends to zero fora=o, 
then the function 


Cy ifs Vaan (18) 


has the value one for u = 0 and decreases to zero as u tends to in- 
finity, having an inflection point for u=an,. If u(w) is of the form 
(11) or of a similar form, then 7(w) is zero for w = 0 and tends to 
1, which value it reaches asymptotically. If, as a very crude ap- 
proximation, we put, with y as a constant, 


T(w) =yw (19) 
with the range 0 < w < 1/y, then 7, and T, become proportional to 


the corresponding average per capita wealth w, and Ww, in the two 
societies. Then we find 


G.= (1l—a)y wm, (20) 


G, Oy We ? (21) 
and 
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G. a Ws 
Tac :: (22) 
G, ib “esau 84 Wy 


Since, as we have said, w. > w, , the ratio G./G, is now even greater 
than the one calculated from equations (6) and (7). If we assume 
that the average per capita wealth is proportional to Y/X , which is 
certainly not true, we find, approximately, 


We a 
RE aes ? (23) 
Wy, TL msnCL 
and 
G, a? 
SS eee (24) 
G. (1 —«q)? 


With this assumption G,/G. could easily be of the order of 10+. Ac- 
tual ratios observed between per capita incomes in western indus- 
trialized countries and the “underdeveloped” countries are of the 
order of 10 or higher. Thus even if G./G, is not as large as 104, it 
may well be of the order of 10°. 

The situation may be complicated further by studying the fol- 
lowing two effects: 


a) The effect of a large G on w. The quantity w is itself an 
increasing function of G. 

b) The effect of what may be called the imitation of past gen- 
erations. The scientific output of the past generation may serve as 
a stimulus toward scientific activity in the present generation. 


Thus we see that two societies, having the same distribution of 
psychological traits and abilities, may differ tremendously in their 
scientific and engineering progress. : 

There are, of course, a number of other factors involved. For 
example, we may consider the effect of the total population N . Con- 
sider two societies both in configuration R, and having the same G. 
If the population N of one society is sufficiently large, so that 
GN >> 1, it will have at any time a sufficient number of scientists 
and inventors. But if the population N’ of the other society is so 
small that GN’ << 1, then there will be only one scientist or inven- 
tor produced in so many years. If the average time 7 between the 
appearance of such an isolated genius is much longer than the aver- 
age life span, then there may be difficulty in transmitting results 
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from one scientist to the succeeding one, even though written records 
may exist. A direct effect of a personal teacher-pupil contact does 
not exist. 

Now let us return to our original two societies I and II. The 
quantities Y, and Y, are determined by the intersection of the F'(w) 
curve (loc. cit.) with the straight line 


k 
ay (25) 
NA 


The constant k is a purely physiological constant and will, therefore, 
be assumed independent of N. But A is partly determined by the 
neurobiophysical parameters, partly by the effectiveness of communi- 
cation between the individuals in a society (Rashevsky, 1950a). To 
imitate someone we must have information as to how that individ- 
ual behaves. Direct information through personal contact, in gen- 
eral, increases with increasing density of population. If S denotes 
the area occupied by the society, then A may be put proportional to 
N/S. Denoting by ¢ a constant of proportionality we see that the 
slope of the straight line (25) is equal to 


kS 
cN2 


(26) 


The smaller the slope, the more stable the configuration of either 
Rk, or k,. But the smaller the slope, the less 1 — a. Hence, for the 
same area, the larger population favors “backwardness” in an “un- 
derdeveloped” society and technical progress in a “developed” one. 
The effect of the area, for the same population, is an opposite one. 


Through a process of learning such as that discussed in chapter 
xxii of loc. cit., a transition from behavior R, to R, may result. It is 
somewhat more difficult to conceive of a transition in the opposite 
direction with a mechanism such as described in chapter xxii of loc. 
cit. Such a transition would result in the stagnation of a hitherto 
developing high technological civilization. Perhaps calamities, affect- 
ing all levels of the society and resulting in a kind of general frus- 
tration, might produce such a result. 

In conclusion it should be emphasized that similar effects may 
be obtained if R, and R, denote some types of behaviors other than 
those discussed above. For example, R, may designate a behavior 
which emphasizes acceptance on faith and development of strong 
dogmatic religious beliefs, while R, consists of a logical analysis of 


N. RASHEVSKY 59 


everything. The quantity ¢ in that case may stand for «, — ¢, of the 
preceding paper (Rashevsky, 1951). 
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Behavior based on acceptance on faith is considered as a result of 
a special conditioning, and is assumed to be mediated through a special 
center. Behavior based on rational critical analysis is mediated through 
a different center. The previously developed theory of two mutually ex- 
clusive behaviors is applied to the situation. Making simple assumptions 
about the development of the two centers with age, it is shown how 
equations can be derived which describe the percentage of individuals 
exhibiting one type of behavior or another as a function of time and of 
different biosociological parameters. 


In chapter xxi of our book (Rashevsky, 1950) a neural mechan- 
ism was suggested which might possibly account for the interplay of 
the following two opposite types of behavior. One type is the tend- 
ency to accept everything uncritically on faith, the other is the tend- 
ency to analyze everything rationally. In the neural mechanism out- 
lined therein each of the two behaviors is mediated through a corre- 
sponding neural center. Each center is characterized by a level of 
spontaneous excitation; e«; for the center which mediates the first 
type of behavior, «, for the center which mediates the second type. 


The difference 
, co) = Ef — ‘Ex (1) 


determines the relative frequency of one behavior or the other ac- 
cording to equations developed previously (Rashevsky, 1950). 

Acceptance on faith may in a sense be considered as a condi- 
tioning process. We become conditioned to certain reactions in given 
situations, and we react in this particular fashion without asking 
ourselves whether the reaction is logically justified. Beliefs and prej- 
udices, including various religious views, are frequently of the na- 
ture of such conditioned reflexes. 

If an individual is strongly conditioned at an early age to a spe- 
cific religious belief, then, even if he later develops his critical rea- 
soning abilities to a high degree, he may have difficulty in abandon- 


61 


62 MATHEMATICAL BIOSOCIOLOGY OF BELIEFS AND PREJUDICES 


ing his early established religious beliefs. Devout Catholics among 
scientists are not rare. Similarly, a man trained in critical reason- 
ing in his own field of investigation may possess racial prejudices 
which were acquired at an early age. 

A certain type of belief or prejudice may be more strongly con- 
ditioned than another. The effect of the center for rational reason- 
ing may be opposed not to a general tendency to accept everything 
on faith, but to a particular situation, accepted on faith by condition- 
ing. Thus for each situation we have an ¢,, which measures the 
tendency to accept this situation arationally. 

It is rather plausible to assume that at birth the individual’s 
center for analytical reasoning practically does not function, so that 
at birth «, = 0. The center for conditioning and for learning certain 
things mechanically and arationally (Rashevsky, 1950) does, how- 
ever, function at a very early age, for the child begins to learn some 
simple reactions very early. We shall, therefore, consider that at 
birth an individual is endowed with a value or ability to learn which 
may be measured by a quantity «,, and that this value « is distrib- 
uted according to some distribution function V,(e). If N is the 
total number of newborn individuals, then 


NV, (&)d €o (2) 


represents the number of newborn individuals who have an «, be- 
tween g& anda +de. 

The quantity e, increases with age, and, for simplicity, we shall 
consider a linear increase. If 4 and +* denote two individual para- 
meters and 7 the age, then we have 


&p = A(r—7"). (3) 


The quantity 4 is, in general, different from individual to indi- 
ca and is distributed according to some distribution function 
Vo(A). 

The same holds about -*, and we should introduce a distribution 
function for that quantity too. For simplicity, we shall consider here 
7 to be the same for all individuals, thus assuming that the individ- 
ual variations in the rate of mental development are due to varia- 
tions of 4 only. The reader will as generalize the a equa- 
tions to include the variability of 7°. 

As regards «;, it may be reasonably assumed to increase with 
-age in proportion to the innatee,. If we denote by a a coefficient 
and by 7) the age at which the conditioning of a certain faith reac- 
tion begins, then for that particular reaction we have 
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Ef = eo (t — 7). (4) 


It must be emphasized that while -* is essentially a biophysical 
constant, 7) is a purely sociological one. 


At a given moment ¢ let the distribution of the quantity ¢ in 
the population be given by a function V(¢,t). From (1), (3), and 
(4) we have 

p= &(1— 7%) —A(r— 7°). (5) 

Let N be the total number of newborn individuals at the time 

t=—0. Of those N individuals the number 
Naas GO) es dh (6) 


will have values of « and /4 in the intervals «=, e& + de, and 
A,4+d21. Atamoment t > 7+ and t > 7* all those individuals will 
have a ¢ given by expression (5) in which we put r= ¢. The num- 
ber of individuals whose ¢ at the time ¢ is between ¢ and ¢ + d¢ is 
equal to the number of individuals whose 4 lies between 


@ &)(t — to) = 
——— a 


b= = 7) 
and 
A—dih, 
where 
aia, (8) 


Introducing (7) and (8) into (6) we find that the number of 
individuals with a given «) and a ¢ between ¢ and ¢ + d ¢ is given by 


Q &o>(t — 70) =p 
t—7 


The total number of individuals at the time ¢ having any value of 
€) and a ¢ in the interval ($, 6 + d¢) is given by 


"ig 0 = 
a ag | vatenve( ot) den (10) 


1 
ie ded. (9) 


aT 


NV; (0) V2 ( 


t—7* eric 


Let N(7)d7 individuals be born in the interval of time 7, rie dr. 
At the time t their ages will be ¢ — +. Hence we see from (10) that 
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at the moment ¢ the number of those individuals with a given ¢ is 
equal to 


: : Gas Erg) = 
ey reek ACD eee (11) 


GG =—& 0 (Se 
, 


If we integrate expression (11) with respect to + between the 
limits — co and t we find the total number of individuals born before 
the moment t and who at the time ¢ have a ¢ between ¢ and¢ +d¢. 
Hence that number is given by 


V(¢,t)d¢ 
(12) 


oO a N (r) h-8)(t — 7. t)i- > 
=d¢ seer AMIS ES deodr. 
t—7—7 t—7—-,7* 


In deriving expression (12) we neglected those individuals who 
were born between ¢ — 7 and t and who have not yet begun to be 
conditioned. If, as is most likely the case, 7* < 7, then, at the mo- 
ment ¢, for all individuals born in the interval (¢ — 7, t — 7”) the 
quantity «, increases, but «; remains zero. For individuals born in 
the interval (t — -*, t) both e«, and «; are zero. For all the individ- 
uals of the first kind ¢; is zero and ¢ is given by 


o——A(b— 7 — 7) or ore (18) 


The number of individuals born between t — rt and ¢ and having at 
the time ¢ a given ¢ is, therefore, equal to 


“N(x = 
ae eri pee dr. (14) 
t-To a at eee ~ 


The total number of individuals who at the time ¢ have a given 


¢, in other words, the expression V(¢, t)d¢, is given by the sum of 
expressions (12) and (14). Hence 
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V(¢,t) 


BG. BNA) @&(t—7—7) —¢ 
a Pate ei ea): Ve Jatesae 

t—r—r t—7;—7* 

0 00 

(15) 

Sie st (7) 

“f =r ee 

t—7r—7* t—7r—7" 

t-To 


The function V(¢é,t) thus contains as a parameter the quan- 


tity +), which is the age at which the conditioning to a certain be- 
lief or prejudice begins. 


We may generalize expressions (12) and (15) by considering 
that the coefficient a varies from individual to individual and that 
its distribution is the same as the distribution of $ at the moment 
7 + t) When the conditioning process begins. This is rather plausible 
because the larger the number of individuals with a positive ¢, the 
stronger will be the effort to condition other individuals toward a 
faith behavior. This stronger effort may be expressed as a larger 
value of a. Instead of (15) we now have: 


Ce le poem (7) Oeht = t= 15) =~ @ 
ven l | a ee) 


a NG) é 
Deg dete da +f alesse 
t—7-— 7" t—7T—T 


t-To 


which is a peculiar type of functional equation for the determina- 
tion of V(¢,t). 

Once the function V(¢,t) is known we also know the percent- 
age of individuals who behave rationalistically and do not accept 
anything on faith. In particular, if the beliefs in question are re- 
ligious beliefs, then the above percentage is also the percentage of 
agnostics in the population. 

Let P.(#) be the percentage of agnostics among all individuals 
with a given $. The function P.(¢) is determined by equations (21) 
and (22) of chapter ii of our book (Rashevsky, 1950). The fraction 
of agnostics p(t) in the whole population at a given time ¢ is 
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p(t) = ie Vie, t)P.le)dg- (17) 


Expression (15) or (16), together with (17), gives us the in- 
cidence of agnostics at different times in terms of such data as the 
distribution of learning ability at birth, the distribution of the con- 
stant 4, which determines the rate of mental development of an in- 
dividual, the birth rates at different times in the past, and, finally, 
the age 7) at which religious conditioning begins. In principle all the 
above functions are determinable experimentally, and thus the equa- 
tions are verifiable. 


It is readily seen that the greater 7), that is, the later the con- 
ditioning begins, the more rapid will be the shift toward a rational 
behavior in the social group. Religious beliefs, which are usually ac- 
quired at a very early age, are most persistent. Certain political be- 
liefs and prejudices, being acquired at a much later age, are less 
stably entrenched in the behavior pattern of society. 


If, as a rough approximation, we consider that an individual 
practically always behaves as an agnostic when his ¢ is below a fixed 
negative value ¢,., then, denoting by U(a) the distribution function 
of a at the time of birth, we find from expression (10) that the num- 
ber N.(7) of individuals who become agnostic at the age z is given by 


N rao ee) oo 
3G lars q fu i Vite.) 


@ &) (7 — 7) p 
x A sce} U ay dist, 


(18) 


Yh ao 


By introducing effects of imitation, as was done in chapter xiii 
of loc. cit., we may study such cases as transition from agnosticism 
to religiousness of individuals who were born and initially raised in 
agnostic surroundings, and for whom initially a = 0. 


The solution of the functional equation (16) is probably diffi- 
cult, but if we use the simpler expression (15), then V(¢é,t) may be 
evaluated by quadratures. We may assume various plausible dis- 
tribution functions for N (x), Valea), anid Vs (2). 


_Instead of equation (8) we may consider the case in which &} 
begins to increase at the age 7* according to 
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d &; 
dt 


ARE SS (19) 
where «x is a constant. Integrated this gives 
A 
é,—— (1—e*#"™) . (20) 
K 


In this case if 7) < 7* and 4 > a e., the quantity ¢ may first be 
positive, then become negative, and then again become positive. Such 
an individual becomes agnostic at a certain age, but later on, at a 
more advanced age, again becomes religious. The smaller x , the more 
advanced the age at which $ again becomes positive; in other words, 
at which the early conditioned beliefs again return. 

As a crude approximation, the following simple treatment of 
some of the above problems may be suggested. Let the total popula- 
tion » increase according to 


dn 

dt 
where # is a constant. Thus, if is the value of n att=—0, 

N= N, e8* , (22) 


Let a fraction a of the newborn individuals eventually become 
agnostic through a mechanism discussed above. And let us further- 
more assume that individuals born of agnostic parents and raised 
as agnostics permanently remain such. If », is the number of ag- 
nostics, then 


AU aren Fy eer (23) 
dt 


Equations (22) and (23) give: 
dn, 

—=afn, e+ B(1—a)m. (24) 
dt 

Integrated with the initial condition n, — 0 for t = 0, equation 


(24) gives 


Sp ge, (25) 


n 


With 6 ~ 10° year* a a ~ 10%, we find a 6 ~ 10. Hence a 
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completely agnostic behavior will be approached by society in about 
10* years. If a ~ 10-4, the necessary time becomes ~ 10-* years. 

If a fraction y of those born in agnostic families eventually be- 
come religious as a result of imitation, then we obtain a nonlinear 
equation 


dn, 
te Te (26) 
A further elaboration of this approach will also have to take 
into consideration the mortality rate. 


LITERATURE 


Rashevsky, N. 1950. Mathematical Biology of Social Behavior. Chicago: Uni- 
versity of Chicago Press. 


